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For more than 50 years we have known that photosynthetic systems harvest solar energy with 
almost unit quantum efficiency. However, recent experimental evidence of quantum coherence dur- 
ing the excitonic energy transport in photosynthetic organisms challenges our understanding of this 
fundamental biological function. Currently, and despite numerous efforts, the causal connection be- 
tween coherence and efficiency is still a matter of debate. We show, through the study of extensive 
simulations of quantum coherent transport on networks, that three dimensional structures charac- 
terized by centro-symmetric Hamiltonians are statistically more efficient than random arrangements. 
Moreover, we demonstrate that the experimental data available for the electronic Hamiltonians of 
the Fenna-Mathew-Olson (FMO) complex of sulfur bacteria and of the crypophyte PC645 complex 
of marine algae are consistent with this strong correlation of centro-symmetry with quantum effi- 
ciency. These results show that what appears to be geometrically disordered complexes may well 
exhibit a hidden symmetry which enhances the energy transport between chromophores. We are 
confident that our results will motivate research to explore the properties of nearly centro-symmetric 
Hamiltonians in more realistic environments, and to unveil the role of symmetries for quantum ef- 
fects in biology. The unravelling of such symmetries may open novel perspectives and suggest new 
design principles in the development of artificial devices. 



I. INTRODUCTION 

The apparatus used by photosynthetic organisms to 
harvest the sun's energy is both complex and highly effi- 
cient. Although it appears differently in different species, 
some features always persist: Photons are absorbed by 
pigments, usually chlorophyll or carotenoid molecules, 
and transmitted to a reaction center (RC), where the 
primary chemical reactions of photosynthesis occur pQ. 
Photosynthetic structures contain far more chlorophyll 
molecules than RCs, and the former constitute molecu- 
lar networks which pass the energy of an absorbed photon 
to the RC, where it is trapped [THS]. 

In recent years, ultrafast optics and nonlinear spec- 
troscopy experiments provided new insight into exci- 
tonic energy transport in photosynthetic organisms (5J— 
[7]. These experiments report the existence of long co- 
herence times, suggesting that quantum coherence may 
play a fundamental role in the highly efficient transport 
of excitations, and trigger renewed theoretical interest. 

A well established strategy to simulate the quantum 
excitation transport in such molecular complexes is the 
dynamical propagation of the electronic Hamiltonian in 
the presence of the many degrees of freedom of the en- 
vironment [8rH5|. This line of thought may suggest that 
proper tuning of the interaction between system and en- 
vironment is responsible for the long coherence times 
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and the efficient transport [16 . However, the compu- 
tationally obtained coherence times (t c ~ 250 fs) [T3] are 
smaller than the experimentally reported (t c ~ 660 fs) 
5|. Moreover, the complexity of the quantum dynamics 
grows exponentially with the number of strongly coupled 
degrees of freedom [15] . rendering approximations un- 
avoidable, even in a full-fledged computational approach. 
Since, in addition, also the available experimental data 
characterizing the complexes have large uncertainties, 
most computational methods use effective descriptions. 
This contrasts with the fact that the non-trivial, system- 
specific information on complex quantum systems is in- 
scribed into the fluctuations of characteristic observables, 
rather than in their average behavior [T?rt20| . 

On the other hand, simplified models [2"2"H2~4"] that 
mimic the dynamics by those of two coupled two level 
systems interacting with the environment, may help to 
identify some minimal ingredients that must be included 
into a general theory. However, because of their simplic- 
ity, they disregard many of the structural features of the 
light harvesting complexes - possibly those that are most 
relevant for the transport efficiency in realistic systems. 

Therefore, in our present contribution, we attempt to 
incorporate most of the documented properties of the 
underlying electronic Hamiltonians, and yet to outline 
a model as simple as possible. We will see that 3D 
random networks are very likely to perform efficiently 
whenever their Hamiltonians exhibit centro-symmetry. 
On this basis, we will provide strong evidence that the 
transport efficiencies across light harvesting complexes 
indeed substantially depend on the structural properties 
of the apparently disordered macromolecular Hamiltoni- 
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ans. In particular, we will show that, specifically for the 
reported Hamiltonians of the FM08 and PC645 pigment 
complexes, centro-symmetry defines an evolutionary ad- 
vantage. 



II. THE MODEL 

Inspired by the structure of the FMO network - seven 
(FM07) [251 [26] or eight (FM08) [273 chromophores that 
are connected through dipolar interactions - we study 
the simplest possible random model which can grasp its 
essential ingredients: A small random network with N 
sites, where the coherent transport of a single excitation 
is generated by the Hamiltonian 



N 
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Here, a±) and a^} mediate excitations and de-excitations 
of sites j and i, from the local electronic ground state 
to the local excited state, and vice versa. The excita- 

(?) (i) 

tion transfer a+ a_ from site i to site j has a strength 
V id = Vj ti = a/rfj, 
interaction, with 
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consistent with an isotropic dipolar 
and the Tj the position 
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vectors of individual sites. Input and output sites de- 
fine the pole of a sphere of diameter d. The positions 
of the remaining molecular sites are randomly chosen 
within this sphere, what induces a random distribution 
of the remaining Vij. Networks for which one of the 
distances r-jj was smaller than d/100 were discarded in 
the analysis, in order to avoid singular coupling strengths 
[2"T] , The excitation initially is injected at the input site 
| in) = |1), from where it is to be transferred to the out- 
put site |out) = \N). The coupling constant Vi,jv = a/d 3 
between these two sites sets the natural time-scale of the 
dynamics induced by H . The intuition is that the addi- 
tional sites, if properly placed, can mediate a multitude 
of transition amplitudes from input to output, which in- 
terfere constructively upon transmission, and thus ease 
the excitation transfer. A network will be considered ef- 
ficient if the initial excitation is transferred to the output 
site in a time significantly smaller than the Rabi coupling 
time T — ■k/2\Vi i n\ between |in) and |out), with high 
probability. For a quantitative assessment, we define the 
figure of merit [28l ES] 



V = max \{N\e 
te[o,T] 



" iH *|l)| 2 , with T - 0.1 x ir/2\V lj2 \ 
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Indeed, we showed earlier that V fluctuates strongly with 
different realizations of the random molecular network, 
and reaches large values close to unity with finite prob- 
ability. However, while it is therefore natural to ask for 
necessary and/or sufficient conditions on the networks's 
structure such as to guarantee large values of V, no de- 
sign principles were so far identified. This is the purpose 
of our present contribution. 



III. RESULTS 

While near-optimal random networks - in the sense 
of giving rise to large values of V ~ 1 - do not exhibit 
apparent symmetries in their 3D geometry, it was ob- 
served that the individual sites' populations indeed do 
display a near-symmetric structure on the time axis, un- 
der exchange of input and output site, as well as of 
properly defined pairs of intermediate sites (see Fig. 2 
in [29 ). This must be inherited from an exchange sym- 
metry of pairs of two-site coupling matrix elements of 
the underlying Hamiltonian, and evokes an analogy with 
N x N centro- symmetric matrices, which are defined 
by Hi j = -Hjv-i+i.jv-j+i, i-e. , JH = HJ, where J is 
the exchange matrix, Jjj = i^jv-j+i! that exchanges 
site 1 with site N, 2 with N — 1, etc. This type of 
Hamiltonians is known to be tunable towards optimal 
excitation transfer [30H32"] . We now quantify the trans- 




FIG. 1: Correlation between efficiency V and centro- 
symmetry e of 10 s randomly generated 3D networks com- 
posed of seven sites. The color code indicates the logarithm 
of the probability density of a pair (log(e), log('P)). The white 
circle close to the top left corner represents a quasi optimal 
configuration found by a genetic algorithm. 

for efficiency ^ of a given Hamiltonian, and correlate 
this with the Hamiltonian's centro-symmetry, measured 
by e — mins \\H — J^ 1 HJ\\ 1 where || . . . || denotes the 
Hilbert-Schmidt norm, i.e., \\A\\ = VTrAfA [33]. The 
quantity e measures the root-mean square deviation of a 
network Hamiltonian from its image under J, minimized 
over all possible permutations S of the intermediate sites 
2, . . . , N — 1. Figure [T] shows the correlation plot on a 
triple logarithmic scale, for N = 7, and sampling 10 s 
independent random realizations of the network. 

A very strong and, given the variance of V and e over 
several orders of magnitude, unambiguous correlation be- 
tween centro-symmetry (small values of e) and trans- 
fer efficiencies is evident. So is the rather low proba- 
bility to generate a high efficiency molecular conforma- 
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tion by random sampling. On the other hand, the ob- 
served strong statistical correlation implies that an it- 
erative search for highly efficient networks in conforma- 
tion space will with very high probability move towards 
higher degrees of centro-symmetry, and hence acts as an 
evolutionary funnel. This is confirmed by a genetic algo- 
rithm [34 that optimizes the structure of the networks to 
achive maximum efficiency: Starting from random net- 
works, located with high probability close to the global 
maximum of the distribution in Figure [l] i.e. approx. at 
log(e) ~ —0.5, log('P) ~ —3, the algorithm generates 
configurations with almost unit efficiency and very low 
centrosymmetry. An exemplary output configuration is 
indicated by the white circle in Figure [T] While the de- 
viation from the centro-symmetry is reduced by approx- 
imately 1.5 orders of magnitude, the network's efficiency 
increases by about 3 orders of magnitude. 

That the observed correlation between V and e is in 
itself of quantum origin is suggested by our initially for- 
mulated intuition on the role of the intermediate sites 
which mediate the excitation transfer from |in) to |out). 
This is underpinned by inspection of the analogously 
generated correlation plots in Figure [2j though in the 
presence of environment-induced dephasing locally at 
each individual molecular site. The excitation trans- 
fer dynamics is then no longer described by the uni- 
tary evolution alone, but rather by a Lindblad equation 
g(t) = — i [H, g(t)] + Ld cp h({>(t)) , with the Lindblad term 
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and |i) and |j) the N-site electronic state where the ex- 



citation is located at the i , respectively j site. As 
the dephasing rate 7 is tuned from zero over one and 
ten to hundred incoherent events per Rabi coupling time 
T, the distribution in the correlation plot essentially ro- 
tates from a strongly increasing efficiency with increas- 
ing centro-symmetry to an essentially centro-symmetry- 
independent distribution, at the largest dephasing rates. 
However, note that the correlation prevails even for rel- 
atively large dephasing rates of 7 = lOT^ 1 , what corre- 
sponds to one incoherent event within the time interval 
that we chose to define V in 

We now need to validate the relevance of our above 
findings for the actual design of light harvesting com- 
plexes as they abound in nature. Clearly, the strong cor- 
relation between transfer efficiency and centro-symmetry 
of the network as evident from Figs. 1 and 2 was de- 
duced from unrestricted, uniform statistical sampling 
over a random Hamiltonian of the form ([lj. However, ac- 
tual biological functional units are sufficiently well char- 
acterized (by spectroscopic means) to locate them in 
a comparably small sub-volume of conformation space. 
Therefore, we now focus on the trade-off between centro- 
symmetry and transfer efficiency using published data 
of the FM08 protein [35] and the PC645 |Mj light har- 
vesting complexes (of sulfur bacteria and marine algae, 
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FIG. 2: Histogram for the open system dynamics of mir- 
rorsymmetry e vs. transport efficiency V, for different val- 
ues of the dephasing rate in the site basis, 7 = (upper left) 0, 
(right) 1, (lower left) 10, (right) 100, in units of the character- 
istic time scale T = 7r/2|VS n ,out|. The color code indicates the 
probability density of a pair (log(e), log('P)). The histograms 
are generated from 10 5 random networks for 7 > 0, and 10 s 
for 7 = (as in Figure [lj but on linear scale). The corre- 
lation between transfer efficiency and centrosymmetry is less 
pronounced for larger values of 7. 



respectively) as our point of departure. Both these com- 
plexes are represented by a network of eight dipolc- 
coupled two-level systems, and thus comply with our ab- 
stract model from above. 

In order to define an unbiased benchmark in compari- 
son to which we can quantitatively assess the centrosym- 
metry properties of the actual FM08 and the PC645 
Hamiltonian, we use an ensemble of random Hamiltonian 
matrices constructed as follows: the off-diagonal elements 
of these random Hamiltonians are modeled as identically 
and independently distributed random variables, which 
are all chosen from the same Gaussian distribution. The 
mean value and standard deviation of this distribution is 
extracted as the mean and as the standard deviation of 
the non-diagonal elements of the published data for the 
FM08 and PC645 complex, respectively (see the comple- 
mentary material in Sec. |VI B ). In order to avoid centro- 
symmetry to be masked by the biological complexes' on- 
site energies induced by their coupling to (environmen- 
tal) background degrees of freedom, we in addition set all 
diagonal matrix elements equal to their average value. 

The distribution of the centro-symmetry e with respect 
to the input and output sites thus obtained is represented 
by the green histograms, for FM08 and PC645, respec- 
tively, in Figure 3, and defines our unbiased benchmark 
distribution. For the FM08, the input site is the site 
recently added to the structure of the protein, and the 
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output site is usually named "site 3" in the biological 
literature [35]. For the PC645, we identify the DBVc 
and PCBd82 molecules with the input and output sites, 
respectively [55] . 
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FIG. 3: Green: Centro-symmetry distributions for Gaus- 
sian random ensembles with mean and standard deviation 
extracted from the original FM08 (left) and PC645 (right) 
structure data [531 I36| . Red: Centro-symmetry distributions 
for efficient systems obtained by evolutionary optimization of 
the FM08 and PC645 Hamiltonians. The vertical lines indi- 
cate the centro-symmetry of the original FM08 and PC645 
Hamiltonians, respectively. 

The vertical lines in Fig. 3 indicate the centro- 
symmetries of the effective FM08 and PC645 Hamilto- 
nians, respectively, extracted from the published data. 
Since, in both cases, the complexes' centro-symmetries 
arc located in the bulk of the unbiased distribution, they 
do not exhibit any significant symmetry properties. This 
is consistent with the time evolution of the single site 
populations as presented in Figure 4, which clearly show 
very unsatisfactory coherent transfer efficiency, for both 
molecular networks. Indeed, as we have checked, all ele- 
ments of the random ensemble display similarly mediocre 
efficiencies. 
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FIG. 4: Single site populations of the FM08 (left) and 
PC645 (right) light harvesting complexes, with inter-site cou- 
pling matrix elements extracted from published structure data 
[351136] . and on-site energies replaced by their mean. Initially, 
the excitation is located at site 1. Clearly, the transfer effi- 
ciency towards site 8 is mediocre, for FM08 as well as for 
PC645. Site labels 1 and 8 correspond to sites 1 and 3 for 
FM08 ,35 , and to DBVc and PCBd81 for PC645 136], re- 
spectively. 



We now assess how close FM08 and PC645 structures 
are, in conformation space, to more centro-symmetric 
conformations, with higher transfer efficiencies. For this 
purpose, we generate a new ensemble of network Hamil- 
tonians, by repeated execution of N=I0000 iterations of 
a simple genetic algorithm (see Sec. VIA), seeded by the 
published network Hamilonians [351 136] as indicated by 
the vertical lines in Figure 3. This results in the red his- 
tograms in Figure 3, which are clearly shifted towards 
higher centro-symmetries (i.e., smaller values of e), by 
more than the original benchmark distributions' widths. 
We verified that all the thus optimized Hamiltonians me- 
diate essentially perfect excitation transfer. Moreover, 
even the average Hamiltonian obtained from an equally 
weighted sum of the optimized ensemble does so! This is 
clearly demonstrated by Figure 5, where the time depen- 
dence of the individual site populations is monitored, as 
generated by this average Hamiltonian. Since different, 
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FIG. 5: Single-site populations for genetically optimized 
FM08 (left) and PC645 (right) light harvesting complexes, 
with inter-site coupling matrix elements given by an equally 
weighted average over 10 3 optimized Hamiltonians generated 
by the genetic algorithm. Site labeling as in Fig. [4] Essen- 
tially optimal transfer efficiencies are achieved in both cases. 

individually optimized Hamiltonians generally do not 
commute, this is a highly nontrivial result which under- 
pins, in particular, that evolutionary optimization seeded 
by the FM08 and PC645 structures converges into essen- 
tially optimal transfer efficiencies, in a statistically robust 
sense. Furthermore, this observation highlights close con- 
nection between centro-symmetry and efficiency, since a 
sum of centro-symmetric Hamiltonians remains centro- 
symmetric. Finally, inspection of the thus obtained av- 
erage optimal Hamiltonian's matrix elements shows that 
most of them can be obtained by weak perturbation of 
those as inferred spectroscopically (see Sec. VI B). 



In contrast to the above results when starting out from 
the FM08 and PC645 structures, the very same genetic 
optimization leaves the distributions as indicated by the 
green histograms in Figure 3 essentially unaffected, with 
correspondingly bad transfer efficiencies, when seeded by 
any of the elements of our random benchmark ensem- 
ble. This markedly distinct behavior of the genetic op- 
timization for FM08 and PC645 on the one hand, and 
of random benchmark Hamiltonians on the other there- 
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fore strongly suggests centro-symmetry as a hidden de- 
sign principle. 

IV. CONCLUSIONS 

Symmetry is a key concept in physics, relevant in fields 
as diverse as High Energy Physics and Condense Matter. 
In our present work, we provide strong evidence to sup- 
port the hypothesis that hidden symmetries in the Hamil- 
tonian of apparently disordered biological structures fos- 
ter quantum capabilities. In particular, we have shown 
that in random artificial networks, but also in the bio- 
logically relevant light harvesting complexes FM08 and 
PC645, (nearly) centro-symmetric Hamiltonians define a 
robust design principle to mediate coherence-induced, ef- 
ficient quantum transport. Note that this hidden symme- 
try property could only be extracted by a statistical anal- 
ysis of suitably constructed matrix ensembles, gauged 
against experimentally inferred average data. This is not 
surprising, since complex quantum systems typically ex- 
hibit very intricate spectral properties, with in general 
only small embedded regions (e.g. on the energy axis 



or with respect to some control parameter (s)) of locally 
well-defined quantum numbers. In contrast to the global 
symmetries of integrable quantum systems, the associ- 
ated (local) symmetries of complex systems are typically 
not given by explicitly defined integrals of motion, and 
this is why they are hard to detect - and yet may manifest 
a dramatic and robust impact on the systems' dynamical 
properties. We trust that our findings may define a good 
starting point for the design of new devices for efficient 
energy transduction [37] , and are confident that the here 
unveiled conspiracy of symmetry and quantum coherence 
in an apparently disordered system may be instrumental 
in other biological systems [3"5] . 
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VI. SUPPLEMENTARY MATERIAL 

A. Genetic Algorithm 

To iteratively optimize the transfer efficiency we proceed as follows: 
Step 1: Read the original Hamiltonian: H*. 

Step 2: Starting from H* , generate M configurations by randomly perturbing all matrix elements H^j of H* , according 
to a Gaussian distribution with standard deviation a, initially set to oo = HZy/lO. 

Step 3: From these configurations (including H*), choose the one with the largest efficiency. This configuration defines 
the new H * . 

Step 4: Repeat steps 2-4 until the maximum number N of iterations is reached. At the nth iteration, perturb the matrix 
element Hij by a Gaussian distribution with standard deviation er„ = ao/n. 

In our optimization of the FM08 and PC645 structures, M — 100 and N = 10000 always achieve convergence. Note 
that the algorithm only optimizes with respect to the transfer efficiency, and not with respect to the cent rosy mmetry. 



B. Hamiltonians 



A comparison of the off-diagonal elements of the average optimized Hamiltonian with the original FM08 and PC645 
Hamiltonians as given in the literature is shown in Figs. [6] and [7j Using the standard notation from the literature, 
the rows in the Hamiltonian of Fig.Jq (FM08) correspond to the sites with labels 8, 1, 2, 4, 5, 6, 7, and 3. Likewise, 
the rows in the Hamiltonian of Fig. (7[(PC645) correspond to sites labeled by DBVc, PCBc82, DBVd, MBVa, MBVb, 
PCBcl58, PCBdl58, and PCBd82. Mean values V and standard deviations Ay which define the random ensembles 
represented by the green histograms in Fig.JSJare derived from the data in Figs. [HJ [7] as: V — 2.6, AV — 34.5 (FM08), 
and V = 21.3, AV = 67.0 (PC645), respectively. 
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FIG. 6: Off-diagonal elements of the average Hamiltonian after optimizing the FM08 structure with a genetic algorithm (left), 
and of the original, spectroscopically inferred FM08 Hamiltonian (right) [35] . 
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FIG. 7: Off-diagonal elements of the average Hamiltonian after optimizing the PC645 structure with a genetic algorithm (left), 
and of the original, spectroscopically inferred PC645 Hamiltonian (right) 36 . 
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